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Abstract 



The two-fractal overlap model of earthquake shows that the contact area distri- 
bution of two fractal surfaces follows power law decay in many cases and this 
agrees with the Guttenberg-Richter power law. Here, we attempt to predict the 
large events (earthquakes) in this model through the overlap time-series analy- 
sis. Taking only the Cantor sets, the overlap sizes (contact areas) are noted when 
one Cantor set moves over the other with uniform velocity. This gives a time se- 
ries containing different overlap sizes. Our numerical study here shows that the 
cumulative overlap size grows almost linearly with time and when the overlap 
sizes are added up to a pre-assigned large event (earthquake) and then reset to 
'zero' level, the corresponding cumulative overlap sizes grows upto some dis- 
crete (quantised) levels. This observation should help to predict the possibility 
of 'large events' in this (overlap) time series. 
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1. Introduction 

The two-fractal overlap model of earthquake [1] is a very recent modeling at- 
tempt. Such models are all based on the observed 'plate tectonics' and the 
fractal nature of the interface between tectonic plates and earth's solid crust. 
The statistics of overlaps between two fractals is not studied much yet, though 
their knowledge is often required in various physical contexts. For example, it 
has been established recently that since the fractured surfaces have got well- 
characterized self-affine properties [2^], the distribution of the elastic ener- 
gies released during the slips between two fractal surfaces (earthquake events) 
may follow the overlap distribution of two self-similar fractal surfaces [1, 5] 
(see also [6]). Chakrabarti and Stinchcombe [1] had shown analytically that 
for regular fractal overlap (Cantor sets and carpets) the contact area distribu- 
tion follows a simple power law decay. 

The two fractal overlap magnitude changes in time as one fractal moves 
over the other. The overlap (magnitude) time series can therefore be studied as 
a model time series of earthquake avalanche dynamics [7, 8]. 

Here, we study the time (t) variation of contact area (overlap) m{t) between 
two well-characterized fractals having the same fractal dimension as one frac- 
tal moves over the other with constant velocity. We have chosen only very sim- 
ple fractals: regular or non-random Cantor sets and random Cantor sets. We 
analyse the time series data of Cantor set overlaps to find the prediction pos- 
sibility of large events (occurrence of large overlaps), we show that the time 
series m{t) obtained by moving one fractal uniformly over the other (with pe- 
riodic boundary condition) has some features which can be utilized to predict 
the 'large events'. This is shown utilizing the discrete or a 'quantum' nature of 
the integrated (cumulative) overlap over time. 

2. The fractal overlap model 

We consider here the overlap sets of two identical fractals, as one shdes over 
the other. We study the overlap distribution P(m.) for the overlap set m defined 
at various (sliding) positions of the two fractals. First, we consider two regular 
cantor sets (at finite generation n). 

(a) (b) 
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Figure 1. (a) A regular Cantor set of dimension In 2/ In 3; only three finite generations are 
shown, (b) The overlap of two identical (regular) Cantor sets, at n = 3, when one slides over 
other; the overlap sets are indicated within the vertical lines, where periodic boundary condition 
has been used. 



Next, we consider the same for two identical random cantor sets (at finite gen- 
erations n). 

(a) (b) 
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Figure 2. (a) A random Cantor set of dimension In 2/ In 3; only three finite generations are 
shown, (b) Overlap of two random Cantor sets (at n = 3; having the same fractal dimension) 
in two different realisations. The overlap sets are indicated within the vertical bars. 

We studied earlier [5] the overlap statistics for regular and random Cantor 
sets, gaskets and percolating clusters [9]. We found a universal scaling be- 
havior of the overlap or contact area (m) distributions P{m) for all types of 
fractal set overlaps mentioned: P'{m') = L'^P{m, L); m' = mL~°', where L 
denotes the finite size of the fractal and the exponent a = 2{df — d); dj being 
the mass dimension of the fractal and d is the embedding dimension. Also the 
overlap distribution P{m), and hence the scaled distribution P'{m'), are seen 
to decay with m or m' following a power law (with exponent value equal to 
the embedding dimension of the fractals) for both regular and random Cantor 
sets and gaskets: P{m) ^ m~/^; (5 = d. However, for the percolating clusters 
[9], the overlap size distribution takes a Gaussian form. 

3. Time series analysis for two Cantor set overlaps 

We consider now the time series obtained by counting the overlaps m{t) as a 
function of time as one Cantor set moves over the other (periodic boundary 
condition is assumed) with uniform velocity. 

Overlap time series data 

The time series are shown in Fig. 3., for finite generations of Cantor sets of 
dimensions In 2/ In 3 and In 4/ In 5 respectively. 
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Figure 3. The time (t) series data of overlap size (m) for regular Cantor sets: (a) of dimension 
In 2/ In 3, at 8th generation: (b) of dimension In 4/ In 5, at 6th generation. The obvious periodic 
repeat of the time series comes from the assumed periodic boundary condition of one of the sets 
(over which the other one slides). 



Cumulative overlap sizes 

Now we calculate the cumulative overlap size Q(t) = mdt and plot that 
against time in Fig. 4. Note, that 'on average' Q{t) is seen to grow linearly 
with time t for regular as well as random Cantor sets. This gives a clue that 
instead of looking at the individual overlaps m{t) series one may look for the 
cumulative quantity. In fact, for the regular Cantor set of dimension In 2/ In 3, 
the overlap m is always 2^, where k is an integer However the cumulative 
Q{t) = Ya=q 2^^' can not be easily expressed as any simple function of t. Still, 
we observe Q{t) ~ ct, where c is dependent on the fractal. This result is even 
more prominant in the case of Cantor sets with df = In 4/ In 5. 




Figure 4. The cumulative overlap Q versus time; for pure cantor sets: (a) of dimension In 2/ In 3 
(at 8th generation) and (b) of dimension In 4/ In 5 (at 6th generation). The dotted line corre- 
sponds to those for two identical but random Cantor sets. In (b) the two lines fall on each 
other. 
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Cumulative overlap quantization 

We first identify the 'large events' occurring at time ti in the m{t) series, where 
m{ti) > M, a pre-assigned number. Then, we look for the cumulative overlap 
size Qi{t) = fl _^ mdt, t < ti, where the successive large events occur at 
times ti-i and U. The behavior of Qi with time is shown in Fig. 5 for regular 
cantor sets with df = In 2/ In 3 at generation n = 8. Similar results are also 
given for Cantor sets with d/ = In 4/ In 5 at generation ra = 6 in Fig. 6. Qi{t) 
is seen to grow almost lonearly in time upto Qi {f i) after which it drops down 
to zero. It appears that there are discrete (quantum) values of Qi{ti). 
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Figure 5. The cumulative overlap size variation with time (for regular Cantor sets of dimension 
In 2/ In 3, at 8th generation), where the cumulative overlap has been reset to value after every 
big event (of overlap size > M where M = 128 (a) and 32 (b) respectively). 




Figure 6. 
In 4/ hi 5 
big event 



Time (t) 



Q 500000 

0) 
N 

W 4OT00 

a 



> 


ID 

.>2 

E, 

o 



(b) 



Time(t) 



1(KK» 12000 i™ 16000 



The cumulative overlap size variation with time (for regular Cantor sets of dimension 
, at 6th generation), where the cumulative overlap has been reset to value after every 
(of overlap size > M where M = 2400 (a) and 2048 (b) respectively). 
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4. Summary and discussion 

If one Cantor set moves uniformly over another, the overlap between the two 
fractals change quasi-randomly with time (see eg., Fig. 3). The overlap size 
distribution was argued [1] and shown [5] to follow power law decay. Here 
we show numerically that if one fixes a magnitude M of the overlap sizes m, 
so that overlaps with m > M are called 'events' (or earthquake), then the 
cumulative overlap Qi{t) = J^._^ mdt, t < ti, (where two successive events 
of m > M occur at times and tj) grows linearly with time up to some 
discrete quanta Qi{ti) = IQo, where Qo is the minimal overlap quantum, 
dependent on M and n (the generation number). Here I is an integer (see Figs. 

5, 6). Although our results here are for regular fractals of finite generation 
n, the observed discretisation of the overlap cumulant Qi with the time limit 
set by n, is a robust feature and can be seen for larger time series for larger 
generation number n. Similar studies for random Cantor set overlap are in 
progress. This model study therefore indicates that one can note the growth 
of the cumulant sesmic response Qi{t), rather than the sesmic event strength 
m{t), and anticipate some big events as the response reaches the discrete levels 
IQq, specific to the series of events. 
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